ADDITIONAL MOON TRACKING 
COMPUTER AND CALCULATOR 

PROGRAMS 



AS-49-17 


DIVISION OF VARIAN 

301 Industrial Way 
San Carlos, California 


4156 



MODIFICATIONS TO HP2000C PROGRAM 


In bulletin AS-49-13, there is a moon tracking program written for the 
HP2000C computer. In that particular program the equations for elevation 
were based on the center of the earth. As a result, there is about a 10 
error in the elevation number. By changing one statement and adding another, 
the error can be reduced to under 0.2°. The following statements show the 
program before and after the changes were made. Also, in statement numbers 
280 and 350 the line reads FOR N=1 to 25. If these statement lines are 
changed to FOR N=1 to 31, then 31 days of data can be requested instead of 
being limited to 25 days. 


BEFORE 


1300 LET H=L6-G 

1310 RFM t CALCULATION OE ELEVATION# E# OF OBJECT 
1320 LET F3=C0S(L5)*C0SCH)*C0S(D1)+SINCD1)*SINCL5> 
1330 LET E2=SQR( l-< E3*E3> > 

1340 LET F=ATNC F3/F2) 

1350 IE E< 0 THEN 18 10 
1360 IF F>I6*R5 THEN 18 10 

1370 REM: CALCULATION OF AZIMUTH# A#OF OBJECT 
1380 LET A2=SI N( D1 >/( C0SCL5) *COSC E) > 

1390 LET A1 = SIN(L5)*SINCD1) + C0S(L5>*C0SCI;1)*C0S(H) 


AFTER 


1300 LET H=L6-G 

1310 REMS CALCULATION OF ELEVATI ON# E# OF OBJECT 
1320 LET F3=COS<L5>*COSCH)*COS(Dl)+SINCDn*SINlL5> 
1330 LET E2=SGRC l-< E3*E3> ) 

1340 LET E=ATN( C E3/E2) -( 1/ ( 6 1 • 33+E2) ) ) 

1345 LET F=ATN( E3/E2) 

1350 IF E< 0 THEN 1810 
1360 IF E>I6*R5 THEN 1810 

1370 REMS CALCULATION OF AZIMUTH# A# OF OBJECT 
1380 LET A2=SIN< D 1 > /( COS < L 5 ) * COS C F ) ) 

1390 LET A2=A2.-( SIN(L5)/C0S(L5) )*( SINC F>/COS( F) ) 
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The following program is a version of Lance Collister's (WB7CCI) original 
moon tracking program written in GE BASIC. This particular version was modi- 
fied by Jay Liebmann, K5JL , to run on a Mits Altai r 8800B. The program re- 

quires about 6700 bytes of memory. For those with a smaller memory, Jay has 
a stripped down version of the program which requires only 4000 bytes. 


500 DIM F(31>*V(31)*Y(31)*Q( 31)*SC31) 

560 P5«2. 0000000000*3. 1415986535 

570 D5* 360*0000000 00 0/P5 

580 R5-P5/360. 0000000000 

582 DKF FNA(X ) * I NT( X* D5* 1 0+ ♦ 5 ) / 1 0 

584 DFF FNBCX )=(X-INT(X) ) *P5 

58 5 PRI NT"WHAT ARE THE STATION CALL LETTERS'*; 

587 INPUT WS 

590 PR I NT "WHAT IS YOUR LATITUDE IN DEGREES* MI NUTES"; 

600 INPUT L5*U5 

610 PR1 NT"WH AT IS YOUR LONGITUDE IN DEGREES* MI NUTES"; 

620 INPUT L6*U6 

630 L5«CL5+U5/60)*R5 

640 L6*(L6+U6/60)*R5 

650 PRI NT"WHAT IS DESIRED PRINTING INCREMENT IN MINUTES"; 

660 INPUT I 

670 PRINT "DO YOU ONLY WANT PRINTOUT WHEN THE MOON IS NEAR THE HORIZON"; 
690 INPUT B$ 

700 IF BS="YES" THEN 730 
710 LET 16=100 
720 GOTO 800 

730 PRI NT"BEL0W WHAT ELEVATION IN DEGREES DO YOU WANT PRINTOUT* 

731 PRINT’TO OCCUR"; 

740 INPUTI6 

750 PRI NT"WHAT ARE THE GMT MONTH* DAY* YEAR DESIRED"; 

760 FOR N»1 TO 31 

770 INPUT FCN)* V<N)*Y(N) 

780 IF F(N) = 0 THEN 860 
78 5 NEXT N 
790 GOTO 760 

800 PRI NT"WHAT ARE THE GMT MONTH* DAY* YEAR* TIME BEGINNING* TIME ENDING 

820 FOR 1 TO 31 

830 INPUT F(N)*V(N)*Y(N)*Q(N)*S(N) 

840 IF F<N)=0 THEN 860 

845 NEXT N 

850 GOTO 820 

860 N5=N- 1 

870 FOR N= 1 TO N5 

880 IF B$= "YES" THEN 900 

890 GOTO 930 

900 E 1=2400 

910 B= 0 

920 GOTO 950 

930 E1 = S<N) 

940 B=Q(N) 

950 M=F(N) 

960 D=V( N) 

970 Y=Y(N) 

980 Yl=Y-( INT( Y/ 100>*100> 

990 PRINT 
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1000 PRINT 

1010 PRI NT"P0SI TI ON OF THE MOON ON" J MJ "/" J DJ "/"J Y 1* "GMT FROM"" "W$ 
1020 PRINT 

1030 PRINT” GMT”*" AZ"*" EL"*" GHA"* " DEC" 

1040 FRINT" "*" --"*” "," " 

1050 PRINT 
1060 I 1 = 2 

1080 I FM>= 3 THEN 1160 

1090 IF INTC (Y- 1853)/4)<1 1 THEN 1120 

1100 C 1 = - 1 

1110 GOTO 1130 

1120 C1 = 0 

1130 J1=365*CY- 1853)*D*30*CM*9)*INT< CM* 10)/ 2) 

1140 J2=INT< CY- 1853)/4)*1*C1 
1150 GOTO 1270 

1160 IF INTC CY- 1852)/4)< 1 1 THEN 1190 
1170 C 1 = - 1 
1180 GOTO 1200 
1190 C1 = 0 

1200 IF M=9 THEN 1240 
1210 IF M* 1 1 THEN 1240 
1220 C2=0 

1230 GOTO 1250 
1240 C2= 1 

1250 J 1 = 365* C Y- 1852>+D*30*CM-3)*INTC CM- 2)/ 2) 

1260 J2=INTC CY- 1852)/4)*C1*C2 

1270 J= J 1 * J2 

1280 T1 = J- 17472- 5 

1290 D9=CB-INTCB/100)* 100)*INTCB/100)*60 

1300 D6=C El -INTC El/ 100)* 100) + INTC El/ 100) *60 

1310 D7«=D9-D6 

1320 D8=D7- I 

1330 IF D7>0 THEN 1350 

1340 GOTO 1380 

1350 IF D8> = 0 THEN 2220 

1360 B*E 1 

1380 T*CB-INTC B/ 1 00) * 1 00) / 1 440* I NTC B/ 100)/ 24 
1390 T5*T1*T 

1400 K1=FNBC .75 12 13* -036601 102*T5) 

1410 K2=FNB< -8 225 13*. 03629 16457* T5) 

1420 K3=»FNB( -995766* - 0027 37778 52*T5) 

1430 K4«*FNBC -97427 1+ - 033863 1922* T5 ) 

1440 K5=FNB( - 0 3 1 2525* - 036748 19 57*T5) 

1450 L8*K 1* -658*R5*SINC 2*K4) *6 - 289*R5* SI N( K2 ) 

1460 L8=L8- 1 .274*R5*SINCK2-2*K4)-. 18 6*R5*SINC K3) 

1470 L8 = L8*-214*R5*SINC 2*K2)- - 1 14*R5*SINC 2*K5) 

1480 L8=L8“ • 059*R5*SINC 2*K2-2*K4>- - 057*R5*SI NC K2*K3- 2*K4) 

149 0 K6=K5+ • 659 3*R5* SI NC 2*K4) *6 - 2303+R5* SI N( K2) - 1 - 27 2*R5*SIN< K2-2*K4) 
1500 L7=5 • 144*P.5*SIN(K6) -- 146*R5*SINCK5-2*K4) 

1520 LETD1 = C0S(L7)*SINCL8 ) * - 3978 2 1*SI NCL7 ) * -9 17 463 

1530 LET D1=ATNCD1 /C SQRC 1 -D1 » 2) ) ) 

1531 Gl = 50*-5*( C Dl)/( -792) )*D5 

1532 G2=80+< (Dl)/( -808) )*D5 

1533 G3=14l-5-C CD1)*C -738)*D5) 

1534 G4= 170-5-C ( D1 )*( -857)*D5) 

1540 A2=C0S(L7)*C0S(L8)/C0SC Dl) 

1550 A 1 = C C0SCL7) *SINCL8 )* .9 17463-SINCL7)*-397821)/C0SCD1) 

1560 A=ATN( A 1/A2) 

1570 GOSUB 1870 
1580 R 1 = A 

1590 L 1= • 0657098 22* T 1 
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1600 L=T*24* l*0027 38 + 6*646055+(L.l-INT(Ll/24)*24> 
1610 L» ( L-I NT( L/24 ) * 24 > 

1630 G=<L/24>*P5-R1 

1640 IF G<P5 THEN 1670 
1650 G=G-P5 
1660 GOTO 1710 
1670 IF G< 0 THEN 1690 
1680 GOTO 1710 
1690 G«G+P5 
1710 H=L6- G 

17 30 E3»C0S<L5>*C0S(H)*C0S< D 1 > + SI NC D 1> *SI NC L5) 

1740 E2-SQRC 1-(E3*E3>) 

17 50 E*ATN< <E3/E8)-( l/( 6 1 • 33*E2> ) ) 

1755 F»ATN(E3/E2) 

1760 IF E< 0 THEN 2178 

1770 IF E> I 6*R5 THEN 2178 

179 0 A2«SIN(Dl)/< C0S<L5)*C0S( F) ) 

1800 A2*A2-( SIN(L5>/C0S(L5> >* < SI NC F) /COSC F> > 

1810 A1«SIN<L5)*SINC D1 ) + C0S(L5)*C0S( DU*C0S<H) 
1820 A1»CSIN(H>*C0S<D1) >/SQRC 1-AH2) 

1830 A»ATN( A1/A2) 

1840 GOSUB 1870 

1850 GOTO 2020 

1870 IF A=» 0 THEN 1890 

1880 GOTO 1930 

1890 IF A2< 0 THEN 1910 

1900 GOTO 2010 

19 10 A=* P 5 / 2 

1920 GOTO 2010 

1930 IF A> 0 THEN 1990 

1940 IF A2< 0 THEN 1970 

1950 A=*P5+A 

1960 GOTO 2010 

1970 A=P5+( A-P5/2) 

1980 GOTO 2010 

1990 IF A2*>0 THEN 2010 

2000 As*A+P5/2 

2010 RETURN 

2020 IF (T-I 1 ) >( 2*1 )/1440 THEN 2040 
2030 GOTO 2050 
2040 PRINT 

2050 IF I NT( B+ • 5> >9 THEN 2080 
2060 S M 
2070 GOTO 2142 

2080 IF INTCB+. 5)>99 THEN 2110 
209 0 SS«" •• 

2100 GOTO 2142 

2110 IF I NTC B+ • 5) >999 THEN 2140 
2120 SS* M *' 

2130 GOTO 2142 
2140 SS« ,M * 

2142 Z 1«FNA( A) 

2144 Z2=FNA(E> 

2146 Z3=FNACG) 

2148 Z4oFNA(D1) 

2150 IF Z 4<0 THEN 2163 

2151 IF Z 3<G 1 THEN 2163 

2152 IF Z 3>G2 THEN 2154 

2153 GOTO 2157 

2154 IF Z 3<G3 THEN 2159 
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2155 IF Z3 >G4 THEN 2163 

2156 GOTO 2161 

2157 Y$*"U" 

2158 GOTO 2170 

2159 YS="V" 

2160 GOTO 2170 

2161 YS=*\J" 

2162 GOTO 2170 

2163 YS«" " 

2170 PRINT S$;STRS(INT(B+.5))^Z 1 > Z 2, Z 3* Z 4; Y $ 
2176 I 1*T 
2178 B=B+ I 

2180 Z*(E-INT( B/ 100>* 100J-60 

2190 IF Z<0 THEN 1290 

2200 E3=INTC B/100)* 100+ 100+Z 

2210 GOTO 1290 

2220 NEXT N 

2230 N=*0 

2240 PRINT 

2260 PRINT 

2270 PRI NT"DO YOU WANT MORE INFORMATION"; 

228 0 INPUT DS 

2290 IF D$="YES" THEN 560 

2300 END 

OK 


The following two scientific calculator programs were contributed 
by Shelby Ennis, W4WNH/8 and William Dayton, WA8BAH. 

The two equations used for these programs are: 


Elevation = sin -1 ( (cos(GHA-LONG)* cosLAT* cosDEC) + (sinLAT • sinDEC) ) 
Azimuth = cos" 1 ( ( cos el?cosLA T ~ C tanLAT ‘ tanEL )) 
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I 
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account your antenna height, or "beam bending". 




















